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ABSTRACT 

The Mott-Hubbard transition in the half-filled Hubbard model is studied 
analytically for the paramagnetic ground state and the classical Neel state. 
The single-particle density of states is obtained by calculating the Green's 
function represented by the infinite continued fraction. The paramagnetic 
solution shows that the Mott-Hubbard transition is signaled by both col- 
lapsing Hubbard bands and appearing 5-function peak at midgap, and the 
transition is second order at the ground state. We also provide specific heats 
in metallic regime to demonstrate that our results recover those of infinite 
dimension. 



PACS numbers: 71.30.+h, 71.10.+X, 75.10.Lp 
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The Hubbard model has been studied first by Hubbard[l], Gutzwiller[2], 
and Kanamori[3] as a model for a narrow band system. It has been recently 
revisited in relation to a possible model for high-Tc superconductors[4,5]. Be- 
sides cuprate superconductors, the Hubbard model may also be applicable to 
describing the metal- insulator transition(MIT) in the materials like ^2-^03 [6] 
in which paramagnetic insulator, paramagnetic metal, antiferromagnetic in- 
sulator, and a new phase antiferromagnetic metal[7] are clearly separated in 
the temperature-pressure and temperature-filling factor(y) phase diagram. 
Transitions between these phases [8], mass enhancement [9] due to strong cor- 
relation and magnetic ordering [9, 10] are common phenomena of the Mott- 
Hubbard systems such as Sri-xLaxTiOslll], Ni{Sei-xSx)2[^'2\ as well as 
V2O3. In addition, an interesting search for a superconducting transition in 
the metallic V2O3 is suggested in close relation to a Mott- Hubbard insulator 
shown in the rare-earth oxide superconductors and the strong electron corre- 
lation shown in the heavy-fermion superconductors. [7] These bring renewed 
interest in the Mott-Hubbard transition(MHT), recently. 

The MHT accompanies both electronic and magnetic transitions. Pure 
V2O3 in a paramagnetic metal phase becomes antiferromagnetic insulator 
through the antiferromagnetic ordering of the vanadium spins at a tempera- 
ture around 150K.[6] Above this temperature, paramagnetic metal-insulator 
transition occurs. This transition temperature, however, can be suppressed 
until very low temperature by applying hydrostatic pressure, adjusting the 
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oxygen stoichiometry, and alloying with titanium sesquioxide. [7-9] There- 
fore, both paramagnetic metal-insulator transition due to strong correlation 
and paramagnetic metal-antiferromagnetic insulator transition due to anti- 
ferromagnetic ordering are interesting. These transitions accompany Fermi 
liquid quasiparticle mass enhancement through the strong electron correla- 
tions, which has been studied by many experiments[8,9,ll]. On the other 
hand, no successful theories exist to explain the experimental results except 
the large-d approach[12-15] as far as we know. 

The beauty of large-d approach lies in the ability to solve the Hubbard 
model, in principle, exactly in a well-defined limit in which spatial correla- 
tions do not play an important role. Unfortunately, some degree of insight 
into the local dynamics contributing to the end-results is concealed as a 
component of the numerical self-consistency that is an integral part of this 
approach. In an effort to compliment this procedure, in this Letter we intro- 
duce an extremely simple approach to the dynamics of the MHT by calculat- 
ing the single-particle density of states(DOS) through a continued fraction 
formalism for the Green's function. The formalism is known as projection op- 
erator technique [16], recurrence relation method[17], and dynamical Lanzcos 
method [5]. The method that we use has been previously applied with great 
success to the response function of the electron gas at a metallic density at 
which usual perturbation technique may not be that successful. [18] Prom the 
study on the electron gas, we found that infinite continued fraction provides 
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analytic expression for the dynamics of the electron system by using minimal 
amount of information about the ground-state correlations. 

Here we study the MHT on a bipartite lattice for the paramagnetic ground 
state and the classical Neel state. The paramagnetic metal-insulator tran- 
sition is realized when temperature is above the Neel temperature. Below 
the Neel temperature, however, the MHT docs not occur at a finite on- 
site Coulomb interaction C/ on a bipartite lattice because of appearing anti- 
ferromagnetic long-range order due to perfect nesting. [19] Even though the 
paramagnetic metal-insulator transition is realized only above the Neel tem- 
perature, the physics of the transition in the artificial paramagnetic ground 
state, for instance, the Brinkman-Rice scenario [20], is also interesting to un- 
derstand the strong correlation effect in the MHT in the Hubbard model. 
This problem, especially in a finite dimension, has been considered as one 
of quite difficult tasks in condensed matter physics, because one must treat 
both spatial correlations and local fluctuations equally rigorously. 

As a lot of insight is gained by reproducing complicated results within 
a very simple theory, in this work we show that a very large fraction of 
the essential physics described in the large-d approach such as the Hubbard 
bands, the quasi-particle peak, and the second order feature of the MHT at 
ground state can be reproduced by a simple continued fraction expansion 
about atomic state. This reproduction may be understood by the arguments 
advanced as part of the large-d approach that a great deal of physics in the 
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MHT is dominated by local dynamics. [12-14] 

As for a finite dimension, however, no first principles calculation has been 
done to obtain the single-particle DOS even for the atomic state as far as 
we know. But our analytic solution keeps the dimensional dependence even 
though approximation becomes cruder as one goes to a lower dimension. 
The effect of spatial correlation which brings more prominent dimensional 
dependence to the DOS will be studied in the future. 

In this Letter, we first show that the classical Neel state which is close to 
the real ground state of the half-filled Hubbard model yields Hubbard gap 
for any [/ > in the single-particle DOS. Thus we recover a well-known 
fact that the MHT occurs at U — for the Hubbard model on a bipartite 
lattice. The paramagnetic solution for the DOS shows the Hubbard feature[l] 
and the Brinkman-Rice feature[20] for the MHT simultaneously in a very 
simple analytic form. We also show that the thermodynamics obtained by 
this DOS reproduces that of large-d, which confirms the similarity of the 
approximations used in the large-d approach and our simple method. 

The single-particle DOS is obtained by calculating the one-particle Green's 
function of the fermion operator at the same site, i.e., (^o|{c]o.(^)) Cjo-}|^o) 
where cij^ and Cj> are fermion creation and annihilation operators with spin 
cr at site j, curly brackets means anticommutator, and |^'o) denotes ground 
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state. The DOS Pa{i^) is given by[l] 



where 



p^iuj) = - lim ImG^- V - te), (1) 
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= ■^O'o{z)\z=iw+e, (2) 

where L is the Liouville operator. The superscript (— ) denotes usual notation 
of the advanced Green's function. The Green's function (2) may be repre- 
sented by an infinite continued fraction [5, 16, 17]. Thus the single-particle 
DOS is written as 

= ^r7 lim Yl R'eao(^) \z=ioj+e, (3) 

TTiV e^0+ J 

where 

«o(^) = , ^ Ai — ^ (4) 

^ - ao H , ' A2 — 

z-aiH ^ 

z-a2+ ' ■ 

where = U)/{U, U), K = (/i., /i/)/ U-i)- The inner product 

is defined by {A, B) = (^o|{^, -B^}|^o) where A and B are operators of the 
Liouville space, B'^ is the adjoint of B. 

The Hubbard model which we study here is written as 

<jl>cr ja 



where < jl > means nearest neighbor sites. To obtain the DOS Pa{i^) where 
a denotes spin-up state, we choose /o = c]^ and calculate a^, and A,^ by using 
the recurrence relation [17] 

f^+i = iLfy - aufu + ^ufu-i (6) 

The key approximation in this process is the Hartree-Fock approximation 
in calculating the ground state correlation function such as (n^- _o.n/^_o.) = 
{nj-a){ni-cr). This corresponds to setting the ground state |^o) as a product 
of atomic states and to neglecting spatial correlations which is valid at large-d 
limit. A similar decoupling scheme has been done on the average containing 
the operator [if, nja\ whose average vanishes, since no steady current exists 
in the ground state. 

First we obtain pcr{oj) for the classical Neel state to study the MHT for 
a state with antiferromagnetic long-range order. We take antiferromagnetic 
spin configuration such that the spin at site j is J,. Then we obtain 

— iU when v is odd, 

= when u is even, (7) 

and 

Ai = qt^, A^ = {q- 1)^2 for z/ > 2. (8) 

To obtain on the other hand, we choose /o = Then we obtain 

tti, = for odd u and iU for even u and have the same A^, as above. Using 
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these, p{uj) = Pa{i^) + P-(t(w) is given by 



Eq. (9) satisfies sum rule exactly and yields Hubbard gap which makes the 
system insulator for any U > 0. We plot Eq. (9) for = and U = It for 
g = 4 in Fig. 1, for instance. 

For the paramagnetic ground state in which one can see the MHT at 
finite U, we obtain the following regular patterns for Aj, and a^, for large-t/ : 

2 

^ \U\t^ + q{q - l)t' 

= — {l + 0(-)}^a, 
A2.+2 = 4gi2{l + 0(^)} = 6, 

where u > I and ao — ai — a2 — iU/2 and — ij{l+0{t'^/U'^)} for u > 3. 
Only even powers oit/U appear in the expansion for /S.^ and a for 1/ > 3. 
A further approximation (q" — 1) ~ ? which is valid at higher dimension has 
also been made. 

By neglecting terms of (9(t^/f/^), the infinite continued fraction (4) can 
be easily calculated, and the result is 



aoU) = 



5 + ^ 



^2 I A2 

^ ^ 26 



{h-a)-z^ + Mz' + a - 6)2 + Ahz'' + Ai 



(10) 



where z = z — ao- We choose positive sign in front of the square root instead 
of negative sign to satisfy the boundary condition ao{t — 0) — 1, where ao{t) 
is the inverse Laplace transform of ao{z). 

If one sets the chemical potential at /i — ^, Eq. (3) gives the DOS as 
follows: 



2b-K\uj\ * 



f (6-a) + Ai + (f -l)a;2 



A 2 



for a>b (11) 



i; + (i-f) 



+ 



207r|a)| 



-2 — ^ — - for a<b, (12) 

t(6-a) + Ai + (f -1)^2] 

where W = {a;^ — {^/a — Vb)^}{{^/a + \/6)^ — cu^}. From the single-particle 
DOS, one can see that MIT occurs at a = 6, i.e., Uc — 4:y/qt and the Fermi 
liquid quasi-particle mode appears at midgap in metallic phase, i.e., U < Uc 
and the Hubbard bands have band width i-^/qt in insulator(a > b) and U 
in metallic phase(a < b). The band gap is \U — 4^i|. To draw Pa{uj) in 
the large-d hmit, we scale t as t = t^/^/2q and set — 1. We plot three 
typical cases, i.e., (a) U = 4(two bands), (b) U ^ Uc ^ 2^2 (single band), 
(c) U — 2.4(two bands with 5-function peak at a; = 0) in Fig. 2. 

The 5-function quasi-particle peak at zero frequency is of course an ar- 
tifact, and a more realistic calculation would lead to a dispersive band. In 
fact, a similar calculation to present one but expanded for small-C/ shows the 
dispersive band, which will be submitted in a separated paper. Nevertheless, 



the weight of the 5-function peak gives us a measure of the quasi-particle 
mass renormahzation Zp — m/m*. The quasi-particle mass diverges as U 
approaches Uc from metaUic side. The Brinkman-Rice scenario for MIT[20] 
can be seen here. 

It is interesting to compare our results to those of large-d. Our DOS 
shows that the MHT in paramagnetic ground state is second order, which 
coincides with the recent results[15] of infinite dimension. The dispersive 
band shown in the large-d results may be reproduced by applying the Lut- 
tinger theorem[21] which fixes the height of the quasi-particle peak to the 
value of the DOS at Fermi energy of the noninteracting system. For further 
comparison with the large-d approach, we extract a part of thermodynamics 
from the DOS (12). 

The extension of our formalism to the finite temperature is rather simple, 
since the temperature-dependence only enters in the correlation function such 
as ijij-^ni-a) which vanishes within Hartree-Fock approximation. Therefore, 
the temperature-dependence does not appear in the DOS in atomic limit. 
Thus the internal energy has temperature- dependence only in the Fermi 
distribution function /(a;). In the metallic regime, the internal energy is 
given by 



where Pmc(<^) is the incoherent part of Eq. (12) and pqp{oj) is the quasi- 
particle DOS. The latter is the coherent part of Eq. (12) divided by the 




(13) 
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quasi-particle mass renormalization Zp, which can be written as 

The parameter rj is determined by the condition Pqp{0) = N{0)/Zf, where 
N(0) = V2/7r[22]. The specific heats calculated from Eq. (13) for C/'s which 
are a little bit less than Uc are shown in Fig. 3. One can see a peak at 
low-temperature regime which is given by the Fermi liquid quasi-particle. 

One of interesting features which looks peculiar among the results of 
largc-d approach is shown in the specific heat at metallic regime(C/ < Uc)- 
The specific heats for different C/'s show a crossing at a certain temperature, 
T 6.5. This phenomenon has also been shown in the experiments for 
^He[23]. But the physical origin of it is not clear yet. These phenomena 
are the characteristics of large-d result [8]. Other thermodynamic quantities 
obtained from Eqs. (11) and (12) will be appeared in the forthcoming paper. 

As final remarks, we mention about our approximation. The Hartree- 
Fock approximation in calculating the static correlation functions is the only 
approximation when we obtain the DOS for the classical Neel state. For 
the paramagnetic solution, however, we neglect 0{t'^/U'^) terms compared 
with unity in obtaining /S.^, and a^, for i/ > 3 in addition to the Hartree-Fock 
approximation. The neglected terms may shifts the positions of poles a little 
bit when they are included. We imagine that this effect may not change the 
over-all physics contained in the DOS. Even under these approximation, our 
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DOS satisfies sum rule exactly, i.e. pa{u!)du! = 1, for [/'s given in Fig. 2. 

The underlying reason for the remarkable agreement between the large- 
d and the simple continued fraction method must lie in the fact that the 
physics of the large-d limit is dominated by the atomic limit. The continued 
fraction expansion used here may provide the possibility to treat the effects 
of quantum fluctuation in the classical Neel state and spatial correlation as 
well as doping. These will be the subjects of forthcoming paper. 
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Figure Captions 

Fig. 1 : The single-particle density of states p{uj) for the classical Neel state. 

The dotted line is for [/ = and the solid hne for U — It. 

Fig. 2 : The single-particle density of states pa{oj) of the paramagnetic 

phase for large-d hmit. (a) U = 4i*, (b) [/ = [/^ = 2%/2U, (c) U = 2AU. 

The number (0.11) in (c) represents the weight of (5-function. 

Fig. 3 : The specific heats for U < Uc — 2v^. The sohd, dashed, and dotted 

hnes are for U—2.6, 2.4, and 2.2, respectively. We set = 1. 
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